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Many weighted scale-free networks are known to have a power-law correlation between strength 
and degree of nodes, which, however, has not been well explained. We investigate the dynamic 
behavior of resource/traffic flow on scale-free networks. The dynamical system will evolve into a 
kinetic equilibrium state, where the strength, defined by the amount of resource or traffic load, is 
correlated with the degree in a power-law form with tunable exponent. The analytical results agree 
well with simulations. 
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I. INTRODUCTION 

A very interesting empirical phenomenon in the study 
of weighted networks is the power-law correlation be- 
tween strength s and degree k of nodes s ^ 
Very recently, Wang et al have proposed a mutual selec- 
tion model to explain the origin of this power-law corre- 
lation 0. This model can provide a partial explanation 
for social weighted networks, that is, although the gen- 
eral people want to make friend with powerful men, these 
powerful persons may not wish to be friendly to them. 
However, this model can not explain the origin of power- 
law strength-degree correlation in weighted technological 
networks. 

In many cases, the concepts of edge-weight and node- 
strength are associated with network dynamics. For ex- 
ample, the weight in communication networks is often 
defined by the load along with the edge @, and the 
strength in epidemic contact networks is defined by the 
individual infectivity On the one hand, although 
the weight/strength distribution may evolve into a sta- 
ble form, the individual value is being changed with time 
by the dynamical process upon network. On the other 
hand, the weight/strength distribution will greatly affect 
the corresponding dynamic behaviors, such as the epi- 
demic spreading and synchronization |a. 1^ ITol . 

Inspired by the interplay of weight and network dy- 
namics, Barrat, Barthelemy, and Vespignani proposed an 
evolution model (BBV model for short) for weighted net- 
works 12, 13]. Although this model can naturally repro- 
duce the power-law distribution of degree, edge-weight, 
and node-strength, it fails to obtain the power-law cor- 
relation between strength and degree. In BBV model, 
the dynamics of weight and network structure are as- 
sumed in the same time scale, that is, in each time step, 
the weight distribution and network topology change si- 
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multaneously. Here we argue that the above two time 
scales are far different. Actually, in many real-life sit- 
uations, the individual weight varies momently whereas 
the network topology only slightly changes during a rela- 
tively long period. Similar to the traffic d yna mics based 
on the local routing protocol 0, 0, IT^ Il7j. we inves- 
tigate the dynamic behaviors of resource/traffic flow on 
scale- free networks with given structures, which may give 
some illuminations about the origin of power-law corre- 
lation between strength and degree in weighted scale-free 
networks. 



II. RESOURCE FLOW WITH PREFERENTIAL 
ALLOCATION 

As mentioned above, strength usually represents re- 
sources or substances allocated to each node, such as 
wealth of individuals of financial contact networks [l8l |. 
the number of passengers in airports of world-wide air- 
port networks |l9j |. the throughput of power stations of 
electric power grids [23, and so on. These resources also 
flow constantly in networks: Money shifts from one per- 
son to another by currency, electric power is transmit- 
ted to every city from power plants by several power 
hubs, and passengers travel from one airport to another. 
Further more, resources prefers to flow to larger-degree 
nodes. In transport networks, large nodes imply hubs or 
centers in traffic system. So passengers can get a quick 
arrival to the destinations by choosing larger airports or 
stations. In financial systems, people also like to buy 
stocks of larger companies or deposit more capital in the 
banks with more capital because larger companies and 
banks generally have more power to make profits and 
more capacity to avoid losses. Inspired by the above 
facts, we propose a simple mechanism to describe the 
resource flow with preferential allocation in networks. 

At each time, as shown in Fig. 1, resources in each 
node are divided into several pieces and then flow to its 
neighbors. The amount of each piece is determined by 
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FIG. 1: Resources in node j are divided into several pieces 
and then flow to its neighbors. The thicker lines imply there 
are more resources flowing. It is worth pointing out that, in 
order to give a clearer illustration we do not plot the resource 
flow into node j or out of node i. 



its neighbors' degrees. We can regulate the extent of 
preference by a tunable parameter a. The equations of 
resource flow are 



leN(j) 



(1) 



where Qj—,i(t) is the amount of resources moving from 
node j to i at time t, Sj{t) is the amount of resources 
owned by node j at time ki is the degree of node i and 
N{j) is the set of neighbors of node j. Hi and j are not 
neighboring, then Qj^i = Qi^j = 0. Meanwhile each 
node also gets resources from its neighbors, so at time 
t + 1 ,W i 



jeN(i) ](^N(i) l€NU) 



(2) 



III. KINETIC EQUILIBRIUM STATE 

The Eq. [3 can be expressed in terms of a matrix equa- 
tion, which reads 



/ail 



S{t + l) = AS{t) := 




121 

where the elements of matrix A are given by 
otherwise 



(3) 



(4) 



Since ^ \o-ij\ — 1, Vj, the spectral radius of matrix A 

i=l 

obeys the equality p{A) < 1, according to the Gershgorin 
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FIG. 2: (Color online) The evolution of the strengths of 
node a and b, where nodes a and b are randomly selected 
for observation. The three cases are in different initial states 
which simply satisfy X^i ~ 1- The exponent a — 1. 




FIG. 3: Scatter plots of Si vs h for all the nodes 



disk theorem Here, the spectral radius, p{A), of a 
matrix A, is the largest absolute value of an eigenvalue. 
Further more, since the considered network is symmetry- 
free (That is to say, the network is strongly connected 
thus for any two nodes i and j, these exists at least one 
path from i to j), A^ will converge to a constant matrix 

for infinite k. That is, if given the initial boundary cou- 
rt 

dition to Eq. |31 (e.g. let ^ Si(0) = 1, where n denotes 
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FIG. 4: (Color online) The correlation between degree and 
strength with different a. In the inset, the relation between 9 
and a is given, where the squares come from the simulations 
and the solid line represents the theoretical result 9 = 1 + a. 

the total number of nodes in network), then Si{t) will 
converge in the hmit of infinite t as lim Si{t) — Si for 

t — *oo 

each node i. 

Consequently, Denote <S := (si, S2 ... Sn)"^, one can 
obtain 

S ^ AS. (5) 

That is, for any i, 

Yl (^"^^/ E ^r). (6) 

From Eq. |S1 it is clear that S is just the kinetic equi- 
librium state of the resource flow in our model. Since 
S = lim A'^S{0), S is determined only by matrix A, if 

k — ^oo 

given the initial boundary condition with S{0) satisfy- 

n 

ing ^ Si(0) = 1. Since matrix A is determined by the 

i=0 

topology only, for each node i in the kinetic equilibrium, 
Si = lim Si{t) is completely determined by the network 

t — >oo 

structure. Si denotes the amount of resource eventually 
allocated to node i, thus it is reasonable to define Si as 
the strength of node i. 

IV. POWER-LAW CORRELATION BETWEEN 
STRENGTH AND DEGREE IN SCALE-FREE 
NETWORKS 

The solution of Eq. (6) reads 



FIG. 5: (Color online) The distribution of strength with 
different a. The inset exhibits the relation between 7 and a, 
where the squares come from the simulations and the solid 
line represents the theoretic analysis 7 = (a + /9)/(l + a). 

where A is a normalized factor. 

In principle, this solution gives the analytical relation 
between and ki when X]jeAf(i) ^7 analytically 
obtained from the degree distribution. For uncorrelated 
networks (22j , statistically speaking we have 

s, = Afci+"^P(A;')fc'", (8) 

k' 

where P{k) denotes the probability a randomly selected 
node is of degree k. Since A^^,, P{k')k'°' is a constant 
when given a network structure, one has Sj ^ kl^"', thus 

s{k) - fci+", (9) 

where s(k) denotes the average strength over all the 
nodes with degree k. 

This power-law correlation s{k) ~ k^ where 9 = 1 + a, 
observed in many real weighted networks, can be con- 
sidered as a result of the conjunct effect of the above 
power-law correlation and the scale-free property. Ob- 
viously, if the degree distribution in a weighted network 
obeys the form P{k) ~ k~^ , one can immediately obtain 
the distribution of the strength 

P{s) ^ (10) 

where the power-law exponent 7 = (a -I- /?)/(! + a). 

V. SIMULATIONS 

Recent empirical studies in network science show that 
many real-life networks display the scale-free property 
[23, thus we use scale- free networks as the samples. Since 
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FIG. 6: (Color online) The distributions of degree (left panel) 
and strength (right panel) with different a. The networks are 
generated by the strength-PA mechanism, and those shown 
here are the sampling of size n = 5000. 



the Barabasi- Albert (BA) model [2J| is the mostly stud- 
ied model and lacks structural-biases such as non-zero 
degree-degree correlation, we use BA network with size 
n = 5000 and average degree {k) — 6 for simulations. 
The dynamics start from a completely random distribu- 
tion of resource. As is shown in Fig. 2, we randomly 
pick two nodes a and 6, and record their strengths vs 
time Sa{t) and Sb{t) for three different initial conditions. 
Clearly, the resource owned by each node will reach a 
stable state quickly. And no matter how and where the 
one unit resource flow in, the final state is the same. 

Similar to the mechanism used to judge the weight of 
web by Google-searching (see a recent review paper psj ] 
about the PageRank Algorithm proposed by Google), the 
strength of a node is not only determined by its degree, 
but also by the strengths of its neighbors (see Eq. 7). 
Although statistically s(fc) ^ for uncorrelated net- 

works, the strengths of the nodes with the same degree 
may be far different especially for low-degree nodes as 
exhibited in Fig. 3. 

In succession, we average the strengths of nodes with 
the same degree and plot Fig. 4 to verify our theoretical 
analysis that there is a power-law correlation s ^ k^ be- 
tween strength and degree, with exponent 9 = 1 + a. Fig. 
5 shows that the strength also obeys power-law distri- 
bution, as observed in many real-life scale-free weighted 
networks. And the simulations agree well with analytical 
results. 



VI. CONCLUSION REMARKS 

In this paper, we proposed a model for resource- 
allocation dynamics on scale-free networks, in which the 



FIG. 7: (Color online) The correlation between degree and 
strength with different a. In the inset, the relation between 9 
and a is given, where the squares come from the simulations 
and the solid line represents the theoretical result 6 — 1 + a. 
The networks are generated by the strength-PA mechanism, 
and those shown here are the sampling of size n = 5000. 



system can approach to a kinetic equilibrium state with 
power-law strength-degree correlation. If the resource 
flow is imbiased (i.e. a = 0), similar to the BBV model 
(12.. , the strength will be linearly correlated with de- 
gree as s{k) ~ k. Therefore, the present model sug- 
gests that the power-law correlation between degree and 
strength arises from the mechanism that resources in net- 
works tend to flow to larger nodes rather than smaller 
ones. This preferential flow has been observed in some 
real traffic systems. For example, very recently, we in- 
vestigated the empirical data of Chinese city-airport net- 
work, where each node denote a city, and the edge- weight 
is defined as the number of passengers travelling along 
this edge per week We found that the passenger 
number from one city to its larger-degree neighbor is 
much larger than that from this city to its smaller-degree 
neighbor. In addition, in Chinese city-airport network 
and US airport network .2] , the power-law exponents are 
9 « 1.4 and 9 w 1.5, respectively, which is within the 
range of 9 predicted by the present model. 

The readers should be warned that the analytical solu- 
tion shown in this paper is only valid for static networks 
without any degree-degree correlation. However, we have 
done some further simulations about the cases of growing 
networks (see Appendix A) and correlated networks (see 
Appendix B). The results are quantitatively the same 
with slight difference in quantity. 

Finally, in this model, the resource flow will approach 
to a kinetic equilibrium, which is determined only by the 
topology of the networks, so we can predict the weight of 
a network just from its topology by the equilibrium state. 
Therefore, our proposed mechanism can well apply to 
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FIG. 8: (Color online) The correlation between degree and 
strength with different assortative coefficients r. The param- 
eter a = 0.5 is fixed. The inset shows the numerically fitting 
value of 9 vs assortative coefficients. The networks are gener- 
ated by the generalized BA algorithm [2^12^ of size n = 5000 
and average degree (Ji) = 6. 



estimate the behaviors in many networks. When given 
topology of a traffic network, people can easily predict 
the traffic load in individual nodes and links by using 
this model, so that this model may be helpful to a better 
design of traffic networks. 
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APPENDIX A: THE CASE OF GROWING 
NETWORKS 

Since many real networks, such as WWW and Internet, 
are growing momently. The performance of the present 
resource-allocation flow on growing networks is thus of 
interest. We have implemented the present dynamical 
model on the growing scale-free networks following the 
usual preferential attachment (PA) scheme of Barabasi- 
Albert Since the topological change is independent 
of the dynamics taking place on it, and the relaxation 
time before converging to a kinetic equilibrium state is 



very short (see Fig. 2) , if the network size is large enough 
(like in this paper n ~ 10'^), then the continued growth 
of network has only very slight effect on topology and the 
results is almost the same as those of the ungrowing case 
shown above. 

Furthermore, we investigate the possible interplay 
between the growing mechanism and the resource- 
allocation dynamics. In this case, the initial network 
is a few fully connected connected nodes, and the re- 
source is distributed to each node randomly. Then, the 
present resource-allocation process works following Eq. 
(2), and simultaneously, the network itself grows follow- 
ing a strength-PA mechanism instead of the degree-PA 
mechanism proposed by BA model. That is to say, at 
each time step, one node is added into the network with 
TO edges attaching to the existing nodes with probability 
proportional to their strengths (In a growing BA network, 
the corresponding probability is proportional to their de- 
grees). Clearly, under this scenarios, there exists strong 
interplay between network topology and dynamic. 

When the network becomes sufficient large (n ~ 10"^), 
as shown in Fig. 6, the evolution approaches a stable pro- 
cess with both the degree distribution and strength dis- 
tribution approximately following the power-law forms. 
Furthermore, we report the relationship between strength 
and degree in Fig. 7, which indicate that the power- 
law scaling, s{k) ~ with 6* = 1 + a, also holds even 
for the growing networks with strong interplay with the 
resource-allocation dynamics. 



APPENDIX B: THE CASE OF CORRELATED 
NETWORKS 

Note that, the Eq. (8) is valid under the assumption 
that the underlying network is uncorrelated. However, 
many real-life networks exhibit degree-degree correlation 
in some extent. In this section, we will investigate the 
case of correlated networks. The model used in this sec- 
tion is a generalized BA model [2^ |27| : Starting from 
Too fully connected nodes, then, at each time step, a new 
node is added to the network and m (< toq) previously 
existing nodes are chosen to be connected to it with prob- 
ability 



fco 



(Bl) 



where pi and ki denote the choosing probability and de- 
gree of node i, respectively. By varying the free param- 
eter fco (> —to), one can obtain the scale-free networks 
with different assortative coefficients r (see the Ref. 
for the definition of assortative coefficient). 

The simulation results are shown in Fig. 8, from 
which one can find that the power-law correlations be- 
tween strength and degree in the correlated networks are 
quantitatively the same as that of the uncorrelated net- 
works, however, the power-law exponents, 0, are slightly 
different. Actually, in the positive correlated networks. 
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the large-degree nodes prefer to connect with some other 
large-degree nodes rather than those small-degree nodes, 
thus there may exist a cluster consisting of large-degree 
nodes that can hold the majority of resource. That clus- 
ter makes the large-degree nodes having even more re- 
source than in the uncorrelated case, thus leading to a 



larger 9. In the inset of Fig. 8, one can find that 9 is 
larger in the positive correlated networks, and smaller 
in the negative correlated networks. However, the ana- 
lytical solution have not yet achieved when taking into 
account the degree-degree correlation, which needs a in- 
depth analysis in the future. 
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